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Abstract. In this tutorial, we introduce the basic concepts and mathe- 
matical tools needed for phase-space description of a very common class 
of states, whose phase properties are described by Gaussian Wigner 
functions: the Gaussian states. In particular, we address their manip- 
ulation, evolution and characterization in view of their application to 
quantum information. 
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1 Introduction 

A Gaussian state is a state with Gaussian Wigner functions [T]. In recent years, an 
increasing attention has been devoted to this class of states, as they play a major role 
in quantum information with continuous variables [1113111] • Besides quantum optics, 
where they are generated with current quantum optical technology, Gaussian states 
naturally appear in the description of optomechanical and nanomechanical oscillators, 
gases of cold atoms and ion traps. Furthermore, Gaussian states exhibit extremality 
properties: among the continuous variable states, they tend to be extremal if one 
imposes some constraints on the covariance matrix [5] . 

The better way to deal theoretically with Gaussian states is to use a suitable 
phase-space analysis. In fact, if the Gaussian character is preserved during the dy- 
namics, it is natural to think of the evolution of a Gaussian state as a transformation 
of the covariance matrix and first-moments vector that fully characterize it. As we 
will see, this is the case if we consider the linear and bilinear interactions used in 
quantum optics labs to generate and manipulate quantum states. Furthermore, the 
Gaussian properties of these states may be preserved also during the dissipative evo- 
lution through noisy channels, both Markovian and non-Markovian. 

The main purpose of this tutorial is to introduce the reader to the phase-space 
description of Gaussian states in view of their applications to quantum information. 
After the definition of Gaussian state and its basic properties in relation with the 
covariance matrix and first-moments vector (Sect. H]), we will focus on the unitary 
evolution through suitable symplectic transformations (Sect.jH]). We will illustrate the 
phase-space approach to describe the generation and manipulation of Gaussian states 
with linear and bilinear interaction of modes (Sect. S]). Single- mode Gaussian states 
will be briefly addressed in Sect. jSl while we will discuss two-mode Gaussian states 
in more details in Sect. [51 In the latter case, the concepts of symplectic eigenvalues 
as well as of standard form of the covariance matrix and local symplectic invariants 
will be introduced and applied to calculate the entropies and the mutual information 
of two-mode Gaussian states and to investigate their separability, entanglement and 
Gaussian quantum discord (Sect. (5]) • Dynamics through Markovian noisy channels will 
be addressed in Sect. [71 and Sect.[8lwill consider the effect of a Gaussian measurement 
performed on a multimode Gaussian state. We will review the main results concerning 
the fidelity between Gaussian states in Sect. [SI Sect. [TUl will close the tutorial with 
some concluding remarks. 



2 Basic notation and definition of Gaussian state 

Each mode fc = 1, . . . , n of a system made of n modes or, more in general, n bosons 
is described by the annihilation and creation operators and aj,, respectively, with 
commutation relations [afe,a|] — Ski- The Hilbert space of the whole system, 7i = 
(81^=1 ^k, is the tensor product of the infinite dimensional Fock spaces J-k of the n 
modes, each spanned by the number basis {\m)k\m&i, i-e., by the eigenstates of the 
number operator a\(Xk- If we address non- interacting modes and use the natural units, 

the free Hamiltonian of the system may be simply written as H = X]fe=i(^l'^fe + \) 
and the corresponding position- and momentum-like operators for the A;-th mode are: 

9fc = (flfc + 4) ' and pk = (flfc - 4) ' (1) 

respectively. The commutation relations [qk^Pi] = iSki associated with qk and pi can 
be rewritten in the following compact form, which will turn out to be very useful for 
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the phase-space analysis: 



RkiRi 



i Oki , (2) 



where R = {qi,pi, . . . ,qmPn)^ is a vector of operators and f2ki — [^]ki are the 
elements of the symplectic matrix: 



k=l ^ ^ 

Note that ff = fl = f2^^ . 

We can now introduce the leading element of this tutorial. A n-mode state de- 
scribed by the density matrix g is a, Gaussian state if its characteristic function: 

x[gKA) =Tt g cxp ^-iA^ilkj (4) 

is Gaussian, namely, if can be written in the following form: 

x[g]{A)^cxp{-lA^f2(Tn^A-iA^n{R)} , (5) 

with A = (ai, bi, . . . , a„, b„)^ e M^" and we defined the covariance matrix (CM): 

aki = [cT]ki = \{Rk.Ri+RuRk) - {Rk){Ri) , (6) 

and (O) — Tr[g O] is the expectation value of the operator O. The vector (R) = 
Tr[pi2] is usually referred to as first-moments vector. The uncertainty relations among 
canonical operators impose a constraint on the CM, corresponding to the inequality 
[ilT]: 

(T+'-f2>0, (7) 

that expresses, in a compact form, the positivity of the density matrix g. 
The exponential appearing in Eq. Q is called displacement operator: 

n 

exp {~iA^f2R} - D{A) = D{X) - (g) Dk{Xk) , (8) 

fc=i 

where A = (Ai,...,A„)^ G C", Afe = ^(afe + ib^), and Dk{\k) = exp{Afc4 - 
A^flfc} are single-mode displacement operators acting on the fc-th mode. Displacement 
operator takes its name after its action on the mode operators, namely: 

D\\)akD{\)^ak + Xk. and D\A) RD{A) = R + A . (9) 

In other words, in phase space composed of the couples of conjugate variables qk and 
Pk, k — 1, . . . , n, it displaces a state by an amount A [5]. 

By Fourier transforming the characteristic function Q, we obtain the so-called 
Wigner function of g pillO): 

W[g\{X) = -1- / d^-A ex^{iA-i2X}x[Q\{A) , (10) 
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with X — {xi, yi, . . . , XniVnY ■ Note that: 

being the 2n-dimensional (5-function. From the identity: 



/ 



.^Mexp{-i^-q^.,..-X}. '^-'""-W"^' . ,12) 



where Q is a real, positive-definite symmetric 2n x 2n matrix, it follows that in the 
case of the Gaussian state (O we have: 

tt" •\/det[crJ 

that is still Gaussian. It is possible to show that pure Gaussian states are the only 
pure states with positive Wigner function [mil2| . 

The same Wigner function (jl3p can be calculated also as follows: 

W[e\{X) = -1- / <f-A exp{zA-X}xb](A) , (14) 

where x[q]{A) = x[q\{^-^), namely: 

x[g]{A) = c:i^i^-\A^<TA-iA^{R)^ . (15) 

The equivalence between Eq. (fTU|) and Eq. (fH)) is due to the fact that an equivalent 
definition of Eq. (fTTj) is: 

j ^^"tI exp{zA^X} ^2"(5(2«)(x)_ (16) 
Eq. is a particular case of the more general s-ordered Wigner function pi[T5] : 

W,[^>](X) - -1- / S-A exp {islAp + zA-r2X} x[q]{A). (17) 

If s = 0, then we have the usual Wigner function ([T^ : if s = — 1 or s = 1 we 
obtain the Husimi Q-function or the Glauber-Sudarshan P- function, respectively [5]. 
Furthermore, by using the relation: 

a r-ordered Wigner function can be transformed into a s-ordered one. The s-ordered 
Wigner function is used to define the nonclassical depth T of a quantum state |14) : 

r=\{i--si (19) 

where s is the maximum value for which Wg \q\{X) becomes positive and semidefinite, 
i.e., a probability distribution. One has T = 1 for number states and T = for 
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coherent states. The nonclassical depth can be interpreted as the minimum number 
of thermal photons that has to be added to a quantum state in order to erase all the 
quantum features of state [13] . 

A Gaussian state is fully characterized by its CM and first-moments vector. For 
instance, the purity fi{g) = Tr[gi^] of the Gaussian state depends only on its CM and 
reads: 

^^'^^ ~ 2" V'DitH ' 
where we used the trace rule in the phase space: 



(20) 



Tr 



Oi O2 



(f )" /^^^ d'^X W[d,]iX) W[d2]iX) , (21) 



which follows from the expansion: 



O = 



(f^x w[d]{x) D{x)nD\x) , 



(22) 



where iT = is the parity operator and D{X)nD\X) = D{2X)n 

UD^ {2X), or, equivalently, by using the characteristic function formalism: 



Tr 



Oi O2 



(27r)" 



which follows from: 



O = 



d^'Mx[0i](A)x[02](-A) 



(23) 



(24) 



We recaU also that Tr[i:>(A)] = (27r)"(5(2")(yl) and Tr[i:>(X)] = {2Tr)"d^^"^ (X). Start- 
ing from Eq. we can also obtain the trace form for the Wigner function: 



W[d] = - Tr[6 D{X)nD^X)]. 



(25) 



Note that the identity operator for n modes has a Wigner function given by 
M^[I](X) = TT-", thus, form Eq. 1^ we have Tr[6] = 2"" d'^^X W[d]{X), from 
which follows the normalization of the Wigner function (IT3l) . The Wigner function for- 
malism allows to easily calculate the expectations of symmetrically ordered products 
of field operators [IS], namely: 



h „k 
t ; 



Tr [e[(at)'^af]J = 1 / d'^X W[g]iX) (al)'^ a 

with, as usual, X = {xi,yi, . . . ,x„,2/„)^, ak = -^{xk + iyk), and: 

8''+'' {x dl + y dt)''+'' 



[0 



dx^dyk {k + h)\ 



(26) 



(27) 



x—y—O 



For the sake of completeness, we observe that the expectations in Eq. ([25]) can be also 
obtained starting from the characteristic function expressed in complex notation: 



Tr [g [(at) 



r\h;,ki 1 _ 



(-1)* 



Qh+k 



(28) 



A=0 
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where = Tr[pZ3(A)] and D{\) has been defined in Eq. ([5]): since its derivatives 

in the origin of the complex plane generates symmetrically ordered moments of mode 
operators, the characteristic function is also known as the moment-generating function 
of the signal g. 

In order to become more familiar with the covariance matrix formalism, we con- 
sider the multi-mode state at thermal equilibrium at temperature T described by the 
density matrix u = ®fe=i '^tii(-^fe) with: 



Tr 



a+a.j (1 -H7Vfe)4«'=+i' 



OO 



m— ^ 



m)kk{m\ , 



(29) 
(30) 



where Pk — fi^k/ (kBT), fee being the Boltzmann constant, and Nk ~ {e^'' — 1)^^ is 
the average number of quanta in the fc-th mode with frequency ujk- Its CM cr^ turns 
out to be diagonal and reads: 



cr^ =0crth(A^fe), 



(31) 



k=l 



where <Tth (Nk) = i(l-H27Vfc)l2 is the 2 X 2 CM of the fc-th single-mode thermal state 
with Nk average photons and Im is the m x m identity matrix. Moreover, recalling 
that {n\D{X)\n) = e~2l-^l _L„(|Ap), Ln{z) being Laguerre polynomials, one can easily 
calculate the expression of the characteristic function of the thermal state (PO)) . that 
turns out to be a Gaussian state. Note that if TV^ — ;> Vfc, then cr^ -> ^l2n, that is 
the CM of the vacuum state of n bosons. 

More in general, the 2n x 2n CM S^, of a n-mode Gaussian state g^, v = {1, . . . ,n}, 
can be re-written in the following block form: 





^12 • 




^12 




■ ^2n 






■ cr„ / 



(32) 



where (Tk and dhk are 2x2 real matrices. In particular, CTk corresponds to the CM of 
the state gk = Tr„\{fcj.[p„] and 6hk is related to the (classical or quantum) correlations 
between the modes h ad k: if 6hk — 0, then ghk — T^ii:v\{h.k}[Qv] — gh ^ gk, that is 
the two modes are uncorrelated, and the CM: 



^hk 



CTh Shk 
Kk ^k 



(33) 



of the state ghk reduces to the direct sum of the two single-mode CMs, namely, 

^hk = Ch ® CTk- 



3 Evolution of Gaussian states 

When an evolution preserves the Gaussian character of a state, it can be described 
with suitable transformations of the position- and momentum-like operators or, equiv- 
alently, of R, that preserve the commutation relations ([2|) . These transformations are 
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called symplectic transformations and are the main tool used to describe the kine- 
matics of Gaussian states in the phase space. 

First of all, we recall that the equations of motion of a classical system of n 
particles described by coordinates {gi, . . . , g„} and conjugated momenta {pi, . . . ,pn} 
with Hamiltonian H can be summarized as: 

Rk - ^, (34) 

where R = (qijPi, . . . , and x denotes time derivative and fl is the symplectic 

matrix defined in Eq. Given a transformation of coordinates R ^ R' = FR, one 
has: 

n TT 

R'.^Fksn^tFu—, (35) 

and thus the equations of motions remain invarians if and only if: 

F flF^ = fi , (symplectic condition) (36) 

which characterizes the symplectic transformations and, in turn, describes the canon- 
ical transformations of coordinates. 

Form the quantum mechanical point of view, a mode transformation R' — FR 
leaves the kinematics invariant if it preserves canonical commutation relations ([2]): 
the 2n x 2n matrix F should satisfy the symplectic condition (155)) . It is worth noting 
that if F and G are symplectic transformations, then also F"^ , F^^ — S2F^ S2^ and 
FG are symplectic: the set of the 2n x 2n matrices satisfying the condition Eq. ([55)) 
form the symplectic group Sp(2n,R). 

An important theorem due to J. Williamson [16] guarantees that every CM can 
be diagonalized through a symplectic transformation [T7]. More in detail, if {dk}^^-^ 
is the set of the symplectic eigenvalues of the 2n x 2n CM cr, namely the moduli of 
the eigenvalues {±dk}'k=i of iflcr, where fl is given in Eq. (|3]), then: 

(T = SWS^, (37) 

where W = ^^^^ dk I2 is a n-mode thermal state with Nk ^ dk — average number 
of photon in the fc-th mode [see Eq. pip ], and S is the matrix which performs the 
symplectic diagonalization (as we will see in the following, if cr describes a physical 
state, then dk > 1/2, and, thus, Nk > 0). Now, since the whole set of the symplectic 
transformations is generated by Hamiltonians which are linear and bilinear in the field 
modes |18,19j, the physical statement implied by decomposition (|57)) is that every 
Gaussian state g can be obtained from a thermal state i' by performing the unitary 
transformation Us associated with the symplectic matrix S, which, in turn, can be 
generated by linear and bilinear interactions. Hence, the density matrix corresponding 
to the decomposition (|37p can be written as: 

g^Usi^ul. (38) 

By using the uncertainty relation ([7]), which is invariant under the symplectic 
group Sp(2n,IR), and the decomposition ([57]) . we have: 

SWS^ + ^f2>0 W> - '-f2 ^ dk > Vfc, (39) 

that is the constraint on the CM by the uncertainty relation leads to the constraints 
dk > 1/2 on its symplectic eigenvalues. From Eq. (PO)) it is straightforward to see that 
a Gaussian state is pure if and only if dfe — 1/2. 
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4 Linear and bilinear Hamiltonians 

In order to preserve Gaussian states, a Hamiltonian should be linear or bilinear in the 
fields mode 1 . This kind of Hamiltonian can be experimentally realized by means of 
parametric processes in quantum optical ^20ii21) . optomechanical [221I23| . microme- 
chanical [H] and cold gases pgll^l^Tll^l^ lgU] systems. Though the actual realiza- 
tion of these transformations necessarily involves parametric interactions in nonlinear 
media, their quantum optical implementation is often referred to as quantum infor- 
mation processing with linear optics, according to the linearity of mode evolution. 
The most general Hamiltonian of this kind can be written as: 

n n n 

^ = E di'^ «I + E ala, + Y: gif &lal + h.c, (40) 

k=l k>l=l k,l=l 

and contains three main building blocks, which represent the generators of the corre- 
sponding unitary evolutions to be described in the following subsections. The mode 
transformation imposed by the Hamiltonian (j40p and, thus, the evolution of the vector 
R and of the CM ex writes: 

R~> FR + d, and cr ^ F cr F'^ , (41) 

where d is a real vector and F a symplectic transformations. Remarkably, the con- 
verse is also true, i.e., any symplectic transformation of the form (|4ip is generated 
by a unitary transformation induced by Hamiltonians of the form (|40)) [18 , 19 . In 
this context, it is worth noting that a useful decomposition of a generic symplectic 
transformation F is the following: 

F = O 1^ O' , (Eulcr decomposition) (42) 

where O and O' are orthogonal and symplectic matrices, while D is a positive di- 
agonal matrix. The physical implication of the Euler decomposition (|42l) is that ev- 
ery symplectic transformation may be implemented by means of two passive devices 
(described by the orthogonal matrices O and O') and by single-mode squeezers (de- 
scribed by D) j31j to be addressed in the following subsections. 



4.1 Displacement operator and coherent states 

The first block of the Hamiltonian in Eq. PO)) contains terms of the form H oc g^^^^ a' + 
h.c. and is linear in the field modes. The corresponding unitary transformations are the 
set of displacement operators we used in Sect. [5] to define the characteristic function. 
The comparison between Eqs. (O and (HTI) shows that the CM is left unchanged by 
the displacement operator while the first-moments vector is displaced. 

The displacement operator is strictly connected with coherent states [5]. For a 
single mode, coherent states are defined as the eigenstates of the mode operator, i.e., 
a\a) = a\a), where a € C is a complex number. Using Eq. ([9]), it can be shown that 
coherent states may be defined also as |a) = D{a)\0), that is the unitary evolution 
of the vacuum through the displacement operator. Properties of coherent states, such 
as the overcompleteness and the nonorthogonality directly follow from those of the 
displacement operator. 

Since the CM of coherent states is the same as the vacuum state one [see Eq. (ISTl) 
with Nk — 0], they are minimum uncertainty states, i.e., they fulfill Ineq. ([7|) with 
equality sign and, in addition, with uncertainties that are equal for the position- and 
momentum-like operators (this is directly seen from the CM). 
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4.2 Free evolution and two-mode mixing 

The second block appearing in the Hamiltonian PO]) . i.e., J2k>i=i 9ki ^l^i^ represents 
two diflterent physical processes. 



4.2.1 Phase shift 

The first process refers to the terms proportional to g^'^^'a^dk and describes the free 
evolution of the modes: in most cases these terms can be eliminated by choosing a 
suitable interaction picture. The effect of free evolution is to add an overall phase shift 
that, for single-mode fields, has no physical meaning, but it is of extreme relevance 
in the case of interference phenomena involving different beams of light, such as the 
interferometric scheme used to implement the homodyne detection. The evolution 
operator may be written as U{9) = exp{— 10 aj,a/c} and acts as a phase rotation on 
the field mode Uk, i.e., U^{0) at U{9) ~ e~'^ afc. Hence, the corresponding symplectic 
matrix reads: 

/ rria ft cin \ 

(43) 

and the evolution of the first-moments vector and the single-mode CM follows form 
Eqs. dH]), with d = 0. 




4.2.2 Two-mode mixing 

The second process, involving different mode operators, describes a linear mixing of 
two modes and, in the quantum optics context, the simplest example corresponds to 
a Hamiltonian of the form H cx a^b + b^a, where for the sake of simplicity we consider 
a system of two modes a = ai and b = 0,2- This Hamiltonian describes the action of 
a beam splitter, i.e., the interaction taking place in a linear optical medium such as 
a dielectric plate. The evolution operator can be recast in the form: 

t/(C) =cxp|Ca^S-Cafo^} , (44) 

where the coupling ^ = </> e*^ G C is proportional to the interaction length (time) and 
to the linear susceptibility of the medium. The two-mode mixer is a "passive" device, 
i.e., the total number of quanta in the two modes is a constant of motion. 
The Heisenberg evolutions of modes a and b are given by: 

C/^ (C) aU{C) ^cos(j)a + e^" sin (j) b , (45a) 
(C) bU{C) ^cos(pb- e-''^ sin (/) a , (45b) 

respectively, and the corresponding symplectic matrix S(^ reads [2]: 

(COS012 sin07?.6/ 
— sm(f>TZg cos(/)l2 

where the 2x2 matrix TZg is defined in Eq. P5|) . The first- moments vector and 
two-mode CMs evolve as usual according to Eqs. (|4T|) . only with d — 0. 
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4.3 Single-mode squeezing 

In the particular case of quantum optics, the last block of the Hamiltonian (j40l) de- 
scribes x^"^^ nonlinear interactions in which a photon in the input (pump) is converted 
into two photons, conserving both the energy and the momentum. If the so-called 
phase-matching conditions are arranged in order to emit the two photons into the 
same mode a, we obtain the single-mode squeezing transformations, which, thus, cor- 
respond to Hamiltonians of the form H cx (a^)^ + h.c. Squeezing has been firstly 
introduced for quadrature operators and refers to a phenomenon in which an observ- 
able or a set of observables exhibit a second moment below the corresponding vacuum 
level [32;. 

The single-mode squeezing operator is usually written as: 

5(e)=exp{i[C(atf -Ta']} , (47) 
which corresponds to the following mode evolutions: 

S\OaS{^)^ cosh ra + e''^' sinh r , (48a) 
S'^ (0 a) S{C} = cosh ra^ + e"*'^' sinh r a , (48b) 

with^ = re^"^. By using the mode transformation in Eqs. (|48l) and the definition of the 
quadrature operators ([1]), it is straightforward to find the symplectic transformation 
describing the single-mode squeezing, namely: 

(cos ■0 sin ^ \ 
(49) 
sin i}) — cos ip I 

which allows to calculate the evolution of the first-moments vector and CM according 
to Eqs. (gl]), but with d = 0. 



4.4 Two-mode squeezing 

The two-mode squeezing transformations correspond to Hamiltonians of the form 
H oc a^w + h.c. and describe x^^^ nonlinear interactions introduced in the previ- 
ous subsection but with the two photons emitted in different modes. The evolution 
operator is usually written as: 

^2(0 =exp{ea^&t . (50) 

The corresponding evolutions of the two modes read: 

slid a S2{C) = cosh r a + e^"^ sinhr fet ^ (5la) 
Sl{i) b 82(0 = cosh r 6 - e'"^ sinh r at , (51b) 

where ^ = re^"^ and the symplectic transformation associated with the two-mode 
squeezer is represented by the block matrix: 

/coshrl2 Ri \ 
^2? = „ „ , , (52) 



\ cosh r 1 



2 



where i?^ is defined as in (HH]) and, as usual, the CM <t and the first-moments vector 
of a bipartite input state transform according to Eqs. (|4T|) with d = 0. 
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Fig. 1. (Color online) Plots of the Wigner function of the single-mode Gaussian state 
Q = D{oi) S{^) iy{Nth) 'S'^(C) D'' (a), for different values of Nth, C a: (a) thermal state, (6) 
squeezed thermal state and (c) displaced squeezed thermal state. 



5 Single-mode Gaussian states 

In the case of a single-mode Gaussian state, Eq. ([55]) may be written as [55] : 

g^Dia)S{0'^t^iN)S\ODHc^), (53) 

and the corresponding CM and first-moments vector can be easily calculated by 
applying the phase-space analysis introduced in Sect. H) In particular the CM can 
be calculated applying the squeezing transformation to the CM crth(-/V) of the 
thermal state i/th(-^), namely or = cr th{N) S'^. The explicit expressions of its 
elements are: 

CTfefe = ii^^ cosh(2r) - (-l)''sinh(2r)cos'0 , (fc = l,2) (54a) 
1 + 2A^ 

0-12 = 0-21 = ^ sinh(2r) sini/' , (54b) 

with ^ — re"^^ , while the first-moments vector reads {R) = \/2(5ie[a], 5m[Q;])^. By 
using Eq. ([20)1 we can calculate the purity of the Gaussian state (|53p . that reads 
fj.(g) = (l-f 27V)-i: the purity of a generic single-mode Gaussian state depends only 
on the average number of thermal photons, as one should expect since displacement 
and squeezing are unitary operations and, thus, do not affect the purity of a state. 
The same observation holds when we address the von Neumann entropy: 

Sv{q) — ~Tr[g\ng] . (55) 

For a single-mode Gaussian state we have: 

Sv{g) - / (Vd^) , (56) 

where: 

fix) = {x + ^) In {x+^)-{x- i) In (x - i) , (57) 



and A/dct [cr] — N + ^ — [2fi{g)] ^ corresponds to the only symplectic eigenvalue of 
the 2x2 positive-definite symmetric matrix cr, since, in this case, the two eigenvalues 
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of iflcr are d± = ±Y^det{CTj, as follows from the Williamson's theorem applied to a 
single-mode CM. 

Starting from Eqs. (|53l) and ([5l)) one can obtain the CMs of the coherent state 
(by setting A'^th = C = 0) or of the squeezed vacuum state (with a = ^ = 0). In Fig. [T] 
the Wigner function of the Gaussian state ([55]) is plotted for different values of the 
involved parameters. 



6 Two-mode Gaussian states 

We can identify different classes of two-mode or, more in general, bipartite Gaussian 
states. However, since these states are the simplest scenario where to investigate the 
fundamental issue of entanglement in quantum information, it is useful to introduce 
equivalence classes of Gaussian states with the same amount of entanglement, i.e., 
Gaussian states whose CMs are connected by local symplectic transformations and, 
thus, are locally equivalent. If we write the CM of a two-mode Gaussian state as: 



A C 
B 



(58) 



where A, B and C are 2x2 matrices, then we can define four local symplectic 
invariants, i.e., quantities that are left unchanged by local symplectic transformations: 

h = det[A], h = dct[B], h dct[C], h = det[cr] . (59) 

The CMs of locally equivalent states can be reduced to the following standard or 
normal form [5^135) : 

(a ci 
a C2 
ci 6 
V C2 b 

where the values of a, b, ci, and C2 are determined by the local symplectic invariants 
namely a? = Ii, b"^ = I2, c\Ci = /s and {ab — cf){ab — Cj) = 14- 
The two symplectic eigenvalues of the CM of a generic two-mode Gaussian state 
can be computed in terms of the symplectic invariants 36 : 



(60) 



..^^ ^'-'^^f^ , (61) 

with A{(t) = h + 12 + 2/3 and, in turn, the uncertainty relation (P^)) reduces to: 

d_ > 1/2 . (62) 

Note that for a pure two-mode Gaussian state we have I4 = 1/16 and A{a) — 1/2, 
i.e., a pure Gaussian state has minimum uncertainty. Moreover, bipartite pure states 
necessarily have a symmetric normal form, i.e., a = 5 in Eq. (j60p . as can be seen by 
equating the entropies of the subsystems. 

A relevant subclass of Gaussian states is formed by the two-mode squeezed thermal 
states (for a general parameterization of an arbitrary bipartite Gaussian state, by 
means of a proper symplectic diagonalization, see Ref. |36l ) , i.e., states generated by 
applying the two-mode squeezing operator to a two-mode thermal state, namely: 

g = S2 vti. (^1 ) ® vtx.{N2 ) si {£) . (63) 
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We can calculate the CM of the state (j63l) as cr = I]2^(Ti^S2^, where is the 
symplectic two-mode squeezing matrix (j52p and cr^ is the CM of the thermal state v 
given in Eq. (j3ip with n = 2. In formula, 

1 f AI2 CRA 

Jf?^ being defined in Eq. ([52|) and, if we assume ^ = r G M: 

A EE A{r, Ni,N2) = {1 + Ni+ N2) cosh(2r) + {Ni - N2) , (65a) 
B = B{r, Ni,N2) ^{l + Ni+ N2) cosh(2r) - (A^i - N2) , (65b) 
C EE C(r, A^i, A^2) = (1 + iVi + N2) sinh(2r) . (65c) 

In particular, if A^i = iV2 = we have so-called twin-beam state (TWB) or two-mode 
squeezed vacuum, that plays a leading role in quantum information with continuous 
variable. The first name, TWB, refers to the fact that it shows perfect correlation in 
the photon number, i.e., it is an eigenstate of the photon number difference a^a — b^b, 
which is a constant of motion as the reader can verify. The second name is instead 
connected with a duality under the action of a balanced beam splitter, since one has: 

f/^f e'') S2{0 e^') = S{^e^') ® S{-^e-^') , (66) 

where U{( — j e'^) is the evolution operator of Eq. (|44)) for a balanced beam splitter, 
S2{S.) is the two- mode squeezing operator (|50p . and >S'(^) is the single- mode squeezing 
operator of Eq. (|47p acting on the evolved mode out of the mixer. In other words, 
a TWB entering a balanced beam splitter evolves into a factorized state composed 
of two squeezed vacua with opposite squeezing phases [37] and, viceversa, a TWB 
may be generated by using single-mode squeezers and a linear mixer as in the first 
continuous variable teleportation experiment |38] . 

6.1 Entropies and mutual information 

5.1.1 Von Neumann entropy 

As we have seen in Sect.[3J a consequence of Williamson's theorem is that every Guas- 
sian state can be generated acting on a thermal state with unitary transformations. 
Thus, the von Neumann entropy of a generic Gaussian state reduces to that of the 
thermal state obtained from it by symplectic diagonalization, since unitary operations 
do not affect the entropy of the whole state. In the case of a two-mode Gaussian state 
Qab with CM (T, using Eq. (|56l) and the additivity of von Neumann entropy ([55]) for 
tensor product states, i.e., SviQA Qb) = Sv{qa) + Sv{qb), we obtain: 

Sv{qab) = f{d+) + f{d^) , (67) 

where /(x) has been defined in Eq. ([57]) and d± are the symplectic eigenvalues of cr 
written explicitly in Eq. (I6ip . 

6.1.2 Mutual information and conditional entropies 

Starting from the von Neumann entropies of the state qab and of the two subsystems 
QA = T^t:b[qab] and qb — T^t:a[qab] it is possible to assess how much information 
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about QAB one can obtain by addressing the single parties. This is of course related 
to the correlations between the two modes and can be quantified by means of the 
quantum mutual information or the conditional entropies |39j. 
The quantum mutual information is defined as: 

hi [qab) = Sv [qa) + Sv {qb) - Sv {qab) , (68) 

and can be easily expressed in terms of the symplectic invariants (|59p of cr and its 
symplectic eigenvalues (|6T|) as follows: 

Im{bab) = f (/h) + / (Vh) - f{d+) - f{d.). (69) 

Note that / = Sv{ba) and / {^/h) = Sv{qb)^ since qa and Qh are a single- 

mode Gaussian states (see Sect. [5]). It is also worth noting that, in the case of pure 
states, the entropies Sy {qa) = Sy {qb) correspond to the unique measure of entan- 
glement for pure bipartite states [3U] . 

The conditional entropies are defined as: 

Sa\b{qab) ~ Sv{qab) — Sv{qb) , 

= f{d+) + }{d^)- fiyh) , (70a) 
Sb\a{qab) = Sv{qab) - Sv{,qa) 

= f{d+) + f{d^)~f[^^, (70b) 

and can also assume negative values. If Sa\b{qab) > 0, the conditional entropy gives 
the amount of information that the party A should send to the party B in order to 
allow for the full knowledge of the overall state qab- If Sa\b{qab) < 0, the party 
A does not need to send any information to the other and, in addition, they gain 
—SA\Bi0AB) bits of entanglement, respectively [analogous considerations hold for 
Ss\AiQAB)]- This has been proved for the case of discrete variable quantum systems 
[41) and conjectured |42) for infinite dimensional ones. 



6.2 Separability of Gaussian states 

A bipartite state qab G Ha <8) "Hb is separable if it can be written as a convex 
combination of product states ^3], namely, qab = J2kPkQ^if^ ^ Qh^^ where pk > 0, 
YlikPk = 1, and Q"^'^ e 'H/i, h — A,B. Finding the convex combination of a separable 
state is a challenging task; nevertheless the separability can be revealed with the aid 
of positive but not completely positive maps. In particular, positivity under partial 
transposition (ppt), that is the transposition applied only to a part of a system, has 
been introduced in entanglement theory by A. Peres '44' as a necessary condition for 
separability. In fact, if we apply, for instance, transposition only to elements of the 

first subsystem Aofa. separable state qab, we have — J2kPk{Q^if^)^ Q^if^ ■ Now, 

\ ~ l^k \ ' transposition corresponds to complex conjugation and the 
transposed matrix is a legitimate density matrix itself, being non-negative, self-adjoint 
and with unit trace. Then none of the eigenvalues of g^^ is negative if g is separable. 
The ppt criterion is usually only necessary and entangled states with positive partial 
transposed density matrix are known to exist and are called bound-entangled states 
|45) . R. Simon, however, has proved that for two- mode Gaussian states it represents 
also a sufficient condition for separability [46] . 
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Since complex conjugation corresponds to time reversal of the Schrodinger equa- 
tion, in terms of continuous variables transposition corresponds to a sign change of 
the momentum variables, i.e., a mirror reflection. For a two-mode system described 
by the density matrix qab , partial transposition with respect to system A will be per- 
formed on the phase space through the action of the matrix Aa = Diag(l, — 1) ® I2, 
where the first factor of the direct sum, representing the mirror reflection, refers to 
subsystem A and the second one to subsystem B (partial transposition with respect 
to subsystem B is obtained in a similar way). Hence, the positivity of the partial 
transposed operator leads to the following uncertainty relation: 

^ + ^f2>0, or (T>~^f2A, (71) 

where a — A a cr A a and flA — A a ^ A a- Furthermore, recalling the definition 
(|59p of the four local symplectic invariants, now we have: 

/i=/i, h^h, h = ~h, h = h, (72) 

where Ik are the symplectic invariants referred to S. Thus, in terms of the symplectic 
eigenvalues d± of the partially transposed CM the ppt criterion reduces to: 



d_>l/2. (73) 



with: 



Z\((t)± A/Z\(cr)2 -4/4 

^ — , (74) 



where A{(t) = Ii + I2 ^ 2/3. Here we have shown that the ppt criterion is necessary 
for separability. As for its sufficiency we refer to the original paper 46 . 

An equivalent necessary and sufficient criterion is based on the evaluation of the 
sum of the variances associated with a pair of EPR-like operators, defined on the two 
different subsystems [4 7) . The insight underlying this criterion is that for an entangled 
state it is possible to gain information on one of the subsystems suitably measuring 
the other one. This criterion leads to an inequality that can be expressed in terms of 
elements of the CM expressed in the standard form namely: 



~ai' + ^~\ci\~\c2\-(l' + ^] <0, (75) 



where we introduced the quantities: 7^ = y (& — l/2)/(a — 1/2), a = 2acosh2ri, 

b — 26cosh2r2, ci = 2ci exp (ri +r2), £2 — 2c2exp[— (ri -1-^2)], and ri and r2 are 
suitable squeezing parameters to transform the CM ((60|) into the so-called Duan 
canonical form (see Ref. [17] for details). A separable state, whether Gaussian or not, 
will not satisfy the above inequality. 

It is worth noting that also the negativity of the conditional entropies (|70p is a 
sufficient condition for entanglement (48) . 



6.3 On the quantification of Gaussian entanglement 



For a two-mode state, a quantitative measure of entanglement can be given on the 
observation that the larger is the violation d- < 1/2 the stronger is the entanglement, 
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or more properly, the stronger the resilience of entanglement to noise [milSOllSTllS^ . 
The logarithmic negativity for a two-mode Gaussian state, is given by [53) : 

S(cr) = max|o,-log2(l_| , (76) 

and it is a simple increasing monotone function of the minimum symplectic eigen- 
value (for < (i_ < 1/2). Thus, it represents a good candidate for evaluating 
entanglement in a quantitative way. 

Another convenient and useful way of looking at the entanglement evolution in 
continuous variable systems is by means of the entanglement of formation (EoF), 
which corresponds to the minimal amount of entanglement of any ensemble of pure 
bipartite states realizing the given state |54[|55| . In general the derivation of an ex- 
pression of the EoF for arbitrary states is not a simple task. 

In the case of a symmetric bipartite Gaussian state with CM given by ((60|) with 
a = b, the EoF reads [55] : 

Ef = f{x,n), (77) 

where f{x) is defined in Eq. (|57p . Xm — {<P^ + l/4)/(2d_), and (i_ is the minimum 
symplectic eigenvalue of the partially transposed CM given in Eq. ([6ip. 

For the two- mode squeezed thermal state ([55)) . in which the standard form of the 
CM is obtained from Eq. (15(11) with a>b and ci = — C2 = c > 0, the EoF is still given 
by Eq. ^ but with ^ : 



(a + b){ab -c^ + \)- 2c^det(cr + ^fl) 
" (a -t- 5)2 - 4c2 ■ ^^^^ 

The EoF of other classes of two-mode Gaussian states can be evaluated by following 
the general prescription proposed in Ref. ^56j . 

Quantitative estimation of entanglement can be also obtained by means of entropy 
functionals [57|- In particular, the degree of entanglement of an ideal bipartite system 
can be assessed following the analysis presented in [55] . 



6.4 Gaussian quantum discord 

The correlations of a bipartite quantum system qab , quantified by the mutual infor- 
mation ([68p . can be divided in a quantum part, known as quantum discord, and a 
classical part [59] . The classical correlations are defined as the maximum amount of 
information we can gain on one part of the system by locally measuring the other 
subsystem, and, thus, can be written as a function of the von Neumann entropies of 
the two subsystems as follows |60j: 

CA\B{gAB) = ma,x{Sv{gA) - ^PkSv{g"iB)} ' C^^) 

k 

where the set {-fffe}. Ilk > and J2k^k — I, represents a positive operator-valued 
measure (POVM), g^i^g — TiBigAB I ® nk]/pk is the conditional state of sub- 
system A when the fc-th outcome occurs in a measurement of subsystem B and 
Pk = T^iAB[gAB I <8) Ilk]- The maximum is taken over all the POVMs performable 
on one subsystem. Classical correlations are thus obtained in correspondence of the 
POVM that minimizes the conditional entropy ^kPk'^v{gA\B)' allowing one to ob- 
tain the highest amount of information on the state of system A. As a matter of fact. 
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the above definition is in general non symmetric with respect to the interchange of 
the subsystems. The quantum discord is then defined as: 

T^a\b{qab) = Im{qab) - Ca\b{qab) , (80) 



Im{qab) being the mutual information ((68)) . 

In the particular case of a two-mode Gaussian state, the Gaussian quantum discord 
is evaluated addressing only Gaussian measurements performed on the subsystems 
and can be written as (for conditional Gaussian measurements on Gaussian states, 
see Sect.© [61,62 : 



'Da\b{qab) = Sv^qb) ~ Sv{eAB) + / (/e^) , (81a) 

= ![J^^-SA\B{eAB). (81b) 

where j{x) has been defined in Eq. ([57]) . Sa\b{.Qab^ is the conditional entropy (|70ap 
and, in terms of the symplectic invariants (l59l) , writes [62] : 



TTimin 



i^3i+v/|-ai-4/4)(/2-i/4) 

2(/2-l/4) 



(/l+4/4)(/2 + l/4)/| - ^ 



(82) 



iih+h-ii-Jiiih+ii-iiY-iiihh 
— — otherwise. 



2/; 

In the case of the squeezed thermal state ((63)) . one has (we set ^ = r) 

1 , 2Ni{l + N2) 



~ 2 l-7Vi+iV2 + (l + 7Vi+iV2) cosh(2r) ' ^ ' 

and the explicit expression of the quantum discord (j8ip can be easily evaluated. 

The quantity / ^y^E"™^ j corresponds to the average von Neumann entropy of 

the conditional single-mode Gaussian state in which is left the subsystem A after 
the Gaussian measurement on subsystem B minimizing the conditional entropy in 
Eq. ([75)) (see Ref. [52] for details of the calculation), ^^^^{qab) can be obtained by 
exchanging the roles of the two subsystems. 

It is worth noting that quantum discord can be nonzero even if the state is sepa- 
rable, which indicates that entanglement is not the only source of quantum correla- 
tions. For instance, there are examples of quantum computational algorithms showing 
a speedup with respect to the classical counterparts, even in the absence of entan- 
glement |63ll64j . States with zero discord represent essentially a classical probability 
distribution embedded in a quantum system, while a positive discord, even on sepa- 
rable (mixed) states, is an indicator of quantumness 65, 66,, and may operationally 
be associated with the impossibility of local broadcasting |62| . 

In the case of pure two-mode Gaussian states, since A{cr) = I1+I2 + 2I3 — 1/2 and 

I4 — 1/16, one has SvigAB) = f (^\J^a\ii) ~ ^' ^^"^ hence the Gaussian quantum 
discord ([5T)) reduces to the entropy of entanglement, i.e., 2?A|B(f?AB) = T^B\AiQAB) = 
SvigA) = Sv{gB)- 



7 Gaussian states in noisy channels 



As one may expect, the dissipative dynamics of a Gaussian states in a Gaussian 
environment, or channel, can be reduced to a suitable transformation of its CM and 
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first-moments vector. In this tutorial we focus on Markovian environments, however, 
it is possible to extend the analysis to non-Markovian ones, as described, for instance, 
in Refs. 68, 69^. 

The dynamics of a single-mode quantum state Qt through a (Markovian) noisy 
environment is governed by the following Master equation: 

Qt = ^{{N + l)C[a] + NC[a)] ~ M*V[a] ~ MV[a)]^ gt , (84) 

where C[6]Qt = iOgtO^ - O^Ogt ~ gtO^O and V[d]gt = 2dgt6 ~ OOgt - gtOO are 
Lindblad superoperators, F is the overall damping rate, while g R and M S C 
represent the effective number of photons and the squeezing parameter of the bath, 
respectively [2]. The terms proportional to C[a] and to C[a^] describe losses and linear, 
phase-insensitive, amplification processes, respectively, while the terms proportional 
to 'D\a\ and I?[a^] describe phase dependent fiuctuations. The positivity of the density 
matrix imposes the constraint |Mp < N{N -f 1). At thermal equilibrium, i.e., for 
M = Q, N coincides with the average number of thermal photons in the bath. 

In order to explicitly derive the evolution of the CM and first-moments vector, we 
transform the Master equation ([5^ into the following Fokker-Planck equation for the 
Wigner function M^[gt](X) associated with gt [2]: 

dtW[gt]{X) = ^(d-xX + d^xrT^dx\w[gt]{X), (85) 



2 

where X = (a;,y)^, dx = {dx,dyY ^ and we introduced the diffusion matrix cr 




_ ^ , 5Re[M] 3m[M] \ 

- {\ + N)-MM]l ■ ^''^ 

If the initial state go is a Gaussian state with CM ctq and first-moments vector 
Xq = Tr[goR], respectively, the Wigner function W{X) of the evolved state under 
the action of the Eq. (|86)) is still Gaussian, but with CM and first-moments vector 
given by (see Ref. [2^ for the explicit calculation): 

o-t = e-"cro + (l-e-^*)(Too, and Xt^g-^/^Xo, (87) 

respectively, which show that (Too is the asymptotic CM when the initial state is 
Gaussian, while Xt = Tr[pti?] is damped to zero. 

The extension to two-mode or, more in general, to n-mode states interacting with 
uncorrelated environments, each described by a Master equation of the form ([M]) . is 
straightforward. In this case, if erg and Xq refer to the CM and first-moments vector 
of the initial state n-mode state and /\. , Nk and Mk are the parameter characterizing 
the environment interacting with the k-th mode, then we have: 

(Tt^Gl^^cToGy^ + {l2n-(Gt)(T^, and Xt ^ G^^Xo, (88) 
where Gt = ®h=i e"'^''*l2 and (Too = ®h=i ^h,oo with: 

^({h + Nn)+MMn] 3m[M,] \ 

'^'''^ V ^HMh] (i + iV„)~5Re[M,]j ■ ^ ' 

Starting form at and Xt^ one can easily evaluate the evolution of all the quan- 
tities addressed in the previous sections, such as purity and, for two-mode states, 
the separability thresholds and entropies (the interested reader can find the explicit 
calculations and a thorough analysis, e.g., in Ref. [2] and references therein). 
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8 Gaussian measurements onto a Gaussian state 

In the previous sections we have reviewed how a Gaussian state can be generated and 
characterized. We have also addressed its kinematics and evolution through noisy 
channels. In order to make this tutorial as complete as possible, we now focus on 
conditional Gaussian measurements [7011711 172 .73 . such as homodyne detection and 
double homodyne detection [2^, performed on Gaussian states (for a PhD tutorial 
on the manipulation of Gaussian states at the photon level see, e.g., Ref. [73]). Con- 
ditional measurements are extremely important in quantum information processing, 
since they are at the basis of quantum teleportation and telecloning protocols and 
allow to generate and manipulate new classes of states |75j . Furthermore, a single 
homodyne detector has been recently used to fully characterize a two-mode squeezed 
thermal state [751177]. 

In order to show how a typical calculation involving Gaussian states and operations 
is carried out, we explicitly derive the characteristic function of a conditional Gaussian 
state obtained by performing a Gaussian measurement on one of its n modes. 

Let us consider the following Gaussian characteristic function with zero first- 
moments vector (extension to non-zero first-moments states is straightforward) as- 
sociated with a n-mode state g [for the sake of simplicity we use the characteristic 
function as defined in Eq. (|15p and drop the explicit dependence on the operators]: 

x(A)=exp{-iA^i:.l} , (90) 

where A = {Ai,A2, . . . , A„)^ € R^" is a column vector and S is the 2n x 2n CM. 
Now we assume to perform a Gaussian measurement on one of the modes, that is a 
measurement described by a POVM with Gaussian characteristic function. Without 
lack of generality, we can assume that the measurement involves mode 1, and, thus, 
the corresponding characteristic function may be written as: 

XuiAi) = TT-i exp {-iyl^cTM Ai - iA^x} , (91) 

where <tm and X are the CM the first-moments vector or, more precisely, the outcome 
of the measurement, respectively. For the sake of simplicity, we write A and the CM 
S in the following block form: 

A=(Ai,A2,...,A„)" = (Ai,A)\ and S = (^^t , (92) 

where A e ^ ]r2 _g ^ r2(«-i) ^ j^2(«-i) ^^^^ symmetric, and C G R^ x R2("-i) , 
making evident the mode undergoing the measurement. The conditional characteristic 
function of the system after the measurement with outcome X is: 

= ^xiAi,A)xMi-A^), (93) 

where we used the trace rule (|23]) and p{X) is the probability of the outcome X: 

piX) = / ^'^^,^.7'^"^ xiAi,A) xm(-Ai) (2^)("-i)5(-A) 
J]R2" (27r)" 

exp{-iX^(A + crM)-ix} 
TT^/det[A + ctm] 



(94) 
(95) 
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where S{—A) = nfc=2 ^^^H"^*;) the product of Kronecker deltas in R^. Note that: 
x(Ai, A) xm(-Ai) = TT-i exp {-^iAi,Afa (Ai, A) + iA'^ x} , (96) 

with: 

/A + (TmC\ , . 

^ = s j ■ 

In order to perform the integral (1931) we observe that cr can be rewritten as follows: 

\^ B-cyA + dMy^c J ' v° i2(n-i) y 

(98) 

The matrix B — C^{A + ctm) is the Schur complement of the matrix cr with 
respect to A + ctm- Now, since: 

M(Ai,A)= (Ai + (A + <tm)-'cA,A), (99) 

Eq. ([931) reduces to: 

X'(A) = -i^ exp {-iA^[B - C^(A + ^m)"'C] A - iA^C^(A + <tm)-'x} 

X J^^^cxp{-^{A'f{A + <TM){A')-t{A'fX}, (100) 

= exp |-iA^[B - C'^(A + crM)"^C] A - iA^C^(A + (tm)"^x} , (101) 

where we performed the change of variables A' — Ai + {A + crM)^^CA. The condi- 
tional state x'(A) is a (n — l)-mode Gaussian state with CM B — C'^{A + crM)~^C 
and first-moments vector C'^{A + crMy^X. Analogously, if we carry out the mea- 
surement on the mode n (actually, the last one), we obtain that the conditional 
state is still Gaussian but with CM A — C{B + (Tm)~^C'^ and first-moments vec- 
tor X^(B + a-uy^C, where, now, A e ]R2(n-i) ^ ]r2(«-i)^ B G ^ R2 and 

C e ]R2("-1) X R2. 



9 Fidelity between Gaussian states 



The fidelity is one of the most important figure of merit in quantum information 
and quantifies the similarity between two states gi and £»2- The Uhlmann's fidelity is 
defined as j78j : 

^(^>l,^?2) = {Tr[(y^^72y^)'/']}^ (102) 

and corresponds to the maximal transition probability between all purifications of 
the two states. 

In the case of single-mode Gaussian states with CMs cr^ and first-moments vectors 
Xfc, A; = 1, 2, Eq. (fTU^ leads to i79j: 



exp {-^ (Xi - X2)^((Ti + (T2yHX, - X2)} 

yATs~vs 



(103) 



with A = det[cri + era] and 5 = 4nLi(det[(Tfc] - \). 
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The problem of finding an analytical formula for the fidelity between n-mode 
Gaussian states has been very recently solved in an elegant way |80j . In particular, 
for two-mode Gaussian states one obtains: 

-^(^?l,^?2) -Tr[gi£)2] (104) 
where X ^ 2VA + 2^/B + \ and: 

dci[flcTiflcT2-\U] ^ _ dct[cri + ^i2] det[<T2 + i^] Q 
det[cri + 0-2] ' dct[cri+cr2] ~ 

are symplectic invariants, fl — u}® uj is the symplectic matrix ^ and, using Eq. (|2ip 
or Eq. ([231), we have: 

Tr[gig2] = . (106) 

A/det[<Ti + (T2\ 

The reader can find the full analysis leading to Eq. (jl04p and the extension to n-mode 
Gaussian states in the original paper [80) . 



10 Conclusions 

In this tutorial we have presented the basic tools and results in order to deal with 
Gaussian states in phase space. We have seen how their generation, manipulation and 
propagation through noisy channels can be described by means of suitable transfor- 
mations of the covariance matrix and first-moments vector. Focusing on two-mode 
Gaussian states, we have addressed their characterization by means of mutual infor- 
mation and conditional entropies, their separability and entanglement properties. As a 
pedagogical example, we have explicitly showed how to calculate the conditional state 
of a multimode Gaussian state that has gone through a Gaussian measurement on 
one of its modes. We have also presented the latest results about the fidelity between 
Gaussian states, that is the most important figure of merit in quantum information. 

After the last meeting I had with prof. Federico Casagrande, I promised him we 
would have discussed my research on Gaussian states and my latest results. Unfortu- 
nately, we had not enough time. Of course, these pages cannot substitute an afternoon 
spent with him, his curiosity and his enthusiasm. . . Nevertheless, I believe he would 
have appreciated this tutorial and I hope it could be a useful tool for students and 
scientists interested in quantum optics, the main topic investigated by Federico during 
his activity at the University of Milano. 



I'd like to thank M. G. A. Paris, F. Benatti and P. Marian for useful suggestions and 
discussions. Financial support from the University of Trieste (FRA 2009) and MIUR (FIRB 
RBFR10YQ3H) is acknowledged. 
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